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Abstract 

Current data from neutrino oscillation experiments are in good agreement with 8 = — | and 8 23 = | under the stan¬ 
dard parametrization of the mixing matrix. We define the notion of “constrained maximal CP violation” (CMCPV) 
for predicting these features and study their origin in flavor symmetry. We derive the parametrization-independent 
solution of CMCPV and give a set of equivalent definitions for it. We further present a theorem on how the CMCPV 
can be realized. This theorem takes the advantage of residual symmetries in neutrino and charged lepton mass ma¬ 
trices, and states that, up to a few minor exceptions, (|d|, 8 , 3 ) = (f, |) is generated when those symmetries are real. 
The often considered p-r reflection symmetry, as well as specific discrete subgroups of 0(3), are special cases of our 
theorem. 
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1. Introduction 

While a coherent picture in leptonic mixing has emerged, important measurements are still lacking. In particular, 
the Dirac CP angle 8 and the exact value of the atmospheric neutrino mixing angle 8 23 are of great interest. Whether 
8 2J is maximal or departs sizably from | has important ramifications for flavor symmetry models [il]. The CP 
phase also has model building impact, and the question of whether the lepton sector violates CP has conceptual 
significance in connection to the matter-antimatter asymmetry via leptogenesis J2l- While maximal atmospheric 
mixing is compatible with data since the observation of atmospheric neutrino oscillations, recently first hints towards 
a Dirac CP angle 8 = -1 have arisen from the appearance and disappearance measurements of T2K [3] when 
combined with reactor antineutrino data. Indeed, global fits iSS confirm a mild preference for this particular 
value of CP phase. 

With these in mind, it is tempting to study the origin of such values of 8 and d 23 within theories of flavor 
symmetry. In particular, the so-called p-r reflection symmetry HiiES was often considered in the literature in 
this respect. It transforms the neutrino fields as (v e , v /( , v,.) —> (v*, v*, v*), leading to |5| = f and ft, 3 = | in the 
standard parametrization of the PMNS mixing matrix [ 1l[ 12|]. In our study, we demonstrate that these two features 
arise as the outcome of “Constrained Maximal CP Violation” (CMCPV), which we will establish in a parametrization- 
independent way by maximizing the Jarlskog invariant under a minimal constraint. 

The framework we will discuss is that a flavor symmetry group G is broken such that the neutrino and charged 
lepton mass matrices are invariant under certain subgroups of G. We will propose and prove a general theorem 
revealing that if the residual flavor symmetries are reafl then the CMCPV is generated. There are a few minor 
exceptions to this theorem which we will clarify in Sec. [3] The p-T reflection symmetry is actually a special case of 
this theorem, which can be shown explicitly after a simple basis transformation. We further deduce some corollaries 
from the theorem which are practically useful in understanding and building models for the CMCPV. For instance, 
specific subgroups of 0(3) can generate CMCPV, so do the models with certain groups under which all neutrino fields 
transform as triplets. As an illustration, we will present a simple model to explicitly realize the CMCPV. 
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This paper is organized as follows. In Sec. [3 we will establish our definition of CMCPV in a parametrization- 
independent way, and give a set of equivalent descriptions. Various physical implications (such as leptonic unitarity 
triangles) from CMCPV are further discussed. In Sec.[3| we present our theorem for the origin of CMCPV and derive 
its corollaries which are important for practical model buildings. We will study applications in Sec.|4| and finally we 
conclude in Sec. [3 Some elaborated mathematical proofs are presented in Appendices A and B. 


2. Parametrization-Independent Formulation of Constrained Maximal CP Violation 


What we mean by “constrained maximal CP violation” (CMCPV) is not merely |d[ = | in the standard parametriza- 
tion of the PMNS matrix Cu.m, but both the |d| = | and 0 p3 = |. In general, a parametrization-independent 
definition of the maximal CP violation should be given in terms of Jarlskog invariant J (130, rather than the CP angle 
5 , because 6 is not rephasing invariant. Furthermore, we will clarify shortly that naively maximizing J without 
constraint is already excluded by experimental data. Hence, introducing the new concept of CMCPV is essential for 
studying the viable maximal CP violation. For this purpose, we hist formulate the CMCPV in a parametrization- 
independent form. 

Definition 1 (CMCPV): 

Constrained Maximal CP Violation (CMCPV) is defined as the maximum of the absolute value of Jarlskog invariant 
under the minimal constraint that the absolute values of the elements in the first row of the PMNS matrix are fixed. 


The Definition 1 is parametrization-independent because it does not invoke any explicit form of the PMNS matrix. 
Note that this is a constrained maximization problem. The Jarlskog invariant J 11311 can be regarded as a function 
of the PMNS matrix U. We are looking for the maximal values of the function J(U ), where U is not an arbitrary 
unitary matrix but a constrained one. We impose this constraint on U by requiring the absolute values of its elements 
in the first row, (\U e \\, \U e z\, |f/ e 3 l), be fixed to certain given values. (Actually, fixing any two of them in the first row 
is enough, due to the unitarity condition |t/ e i| 2 + |C/ e 2 | 2 + |t/« 3 l 2 = 1.) The reason that we fix absolute values of the 
elements of U in its first row, rather than any other rows or columns, will become clear shortly (cf. footnote-2). 

This constraint is necessary because without it / would reach its maximal value as allowed by unitarity, J = ^=, 
which is equivalent to all unitarity triangles being equilateral. The corresponding U in this case is just the Wolfenstein 
mixing matrix 0]. 
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which has been excluded by oscillation data, 
invariant is given by 


In the standard parametrization 112ll, the Jarlskog 


J = — sin 6 cos 0i3 sin 2013 sin 2023 sin 2 On . (2) 
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Indeed, if we compute its maximum by d g J = 0 with 0 = {0 p , 0 13 , 0, 3 , 6), we obtain the Wolfenstein mixing, 
(0 13 , 0 [2 , 0, 3 , |d|) = ^arctan-^L, |, This includes the desired values of (0 23 , 6), but gives unrealistic (0 12 , 0 13 ). 

Hence, the Wolfenstein mixing is already excluded by experimental data. To derive acceptable maximization of J , 
we observe that if we fix 0 p and 0 [3 (to their best-fit values for instance) and then maximize /, we still obtain 
(|d|, 0, 3 ) = (f, |). This is in fact consistent with the above parametrization-independent Definition 1 of CMCPV, 
because fixing (0 I2 , 0 13 ) corresponds to fixing the absolute values of the elements in the first row of U under the 
standard parametrization. Note that this is the allowed minimal constraint we could impose on the Jarlskog invariant: 
fixing any other row or column of the PMNS matrix and then maximizing J will not lead to experimentally acceptable 
results^ Hence, the above Definition 1 gives a minimal parametrization-independent definition of viable maximal CP 
violation. 


2 To be explicit, we have directly verified that fixing the second row or the third row of the PMNS matrix will result in | U e j = \U T j \, or 
\U e j\ = |U j \. (j = 1,2,3), respectively. Fixing the first, second or third columns will lead to \IJ = |t/ f3 |, |t/ f] | = \U (3 \, or |t/ f] | = |t/ p l ■ 
(£ = e, \i, r), respectively. All these cases are already excluded by the current neutrino oscillation data. 
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It is worth noting that the Jarlskog invariant can be further cast into a manifestly parametrization-independent 
form Q 
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The above Eg. (l3aT> is a general expression with (a,(5) = e,p, t and (7, /) = 1,2,3, where a + /J and i + j. The 
Eq. (I3b1 > corresponds to the case of (cr, ft) = (e, /j) and (i, j) = (2, 3). It is clear that the Jarlskog invariant can be 
fully determined by any 4 independent matrix elements {|f/ m |, \U a j\, \Upt\, \Upj\) with a + (5 and i 4 j , as shown 
in Eq. (Bal l. According to our above Definition 1 of CMCPV and using Eq. (I3b1 >. we can maximize Jarlskog invariant 
by imposing the extremal conditions with the two elements (\U e \\, \U e i\) of first row fixed, dJ 2 /d\U M i\ 2 = 0 and 
dJ 2 /d\U^\ 2 = 0. From these two equations, we can derive the solution of CMCPV, 


U J = Kf = \{'-\U ej \ 2 ), C/=1.2,3), 


(4) 


which is the maximum of J under the minimal constraint of Definition 1, and is manifestly parametrization-independent. 
The detail of this derivation is presented in Appendi x A| 

When adopting the standard parametrization lfl2ll of PMNS mixing matrix U, we can use Eq. © to immediately 
deduce the explicit realization of CMCPV, (|d|, 0 23 ) = (f, , which is proven in Appendix A Furthermore, we find 

that once we realize (|<5[, f? 23 ) = under the standard parametrization, the PMNS matrix exhibits an interesting 

feature, which we explain as follows. The standard parametrization of the PMNS matrix is expressed as [12], 


U = 




■ Ct^Sr,^ S ,-,6 


Si-tSryi — c„c~s„e 


— Ci-fSry, — 


(5) 


where we have used the common notations ( s ( .., c ( ..) = (sin 6^, cos 6 j .). Under a rephasing 

U' = diag(l, e~ lS , e~ is ) U diag( 1, 1, e is ), 

we obtain 
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For (6, ft, 3 ) = (+f, we find that U' becomes 
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Eq. m explicitly reveals that the first row is real, while the second row and the third row are complex conjugates 
of each other. We call this feature the “row conjugation equality” (RCE). It is easy to see that the reverse is also 
true: holding RCE will result in (|d|, 0 23 ) = (|, | j. Even though we have demonstrated RCE by using the standard 
parametrization ©0, we stress that the RCE form should be independent of parametrizations (up to trivial rephas¬ 
ing). We can use any other parametrizations [ lji]] [ 17] an d maximize Jarlskog invariant under the same constraint as in 
Definition 1. Then, we find that the mixing matrix always exhibits the RCE form, up to a trivial rephasing. In fact, we 
see that any specific RCE form does obey | U \ = \ U T] \ as in our parametrization-independent general solution (|4]i of 
the CMCPV. 

For later usage, let us introduce the following lemma on the RCE. 

Lemma [ 0(3) invariance of RCE ]: 

If a unitary matrix V has the form of RCE, then after a right-handed real transformation V ~> V’ — VR, the matrix 
V' should still have the form of RCE, where R € 0(3) is an arbitrary orthogonal matrix. 
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The proof of this statement is delegated to Appendix B. 1 


Another interesting feature of RCE concerns leptonic unitarity triangles (LUTs) of the PMNS matrix U, in con¬ 
nection to its column orthogonality. 


U* u ■ + U* U ■ + U* U ■ = 0, 

ei ej fit. fij ti T J 


(9) 


where the column indices i + j. We call these unitarity triangles the column triangles. It is evident that if U has the 
form of RCE, then all column triangles should be isosceles triangles because under the RCE the two sides ULt/ ■ 
and U* Tj U Tj have equal length, It/’,(7^1 = \U* i U T] \. 

Unitarity triangles are intrinsically connected to CP violation because all these triangles have the same area, which 
equals half of the absolute value of Jarlskog invariant J. The LUTs are less studied than the unitarity triangles in 
quark sector since measuring the LUTs and thus the leptonic CP violation is much harder. Nevertheless, the LUTs can 
be directly measured in principle via oscillation experiments 0. Furthermore, the LUTs can provide a geometrical 
formulation of the CMCPV. Since we define CMCPV as J reaching its maximum under certain constraints, it also 
means that the area of the LUT reaches its maximum under those constraints. How do these constraints appear 
in our current geometrical picture? The constraint in our Definition 1 is that the first row of U is fixed, namely, 
\U* e \U e ~\, \U* 2 U e3 \ and |t/* 3 t/ el | are fixed, which means that the e-sides of the column triangles are fixed. Hence, 
the Definition 1 is equivalent to saying that each column triangle reaches its maximal area with its e-side fixed. This 
provides a geometrical formulation of the CMCPV. 

Note that for a triangle with its e-side fixed and its perimeter (the sum of the lengths of its three sides) bounded 
from above, its area reaches the maximum if and only if it is an isosceles triangle. This is clear from geometrical 
intuition. In Ref. 0, we proved that a unitarity triangle must always have its perimeter equal or less than 1. This is 
a necessary and sufficient condition for a triangle to be unitarity triangle, and requires the perimeter of each unitarity 
triangle to be bounded from above, which ensures the area of each unitarity triangle to have a maximum. With these, 
we give a geometrical formulation of the CMCPV: it corresponds to the maximal area of the LUT by fixing its e-side, 
and such LUT is an isosceles triangle. 

Finally, we summarize the analysis of this section into the following theorem. 

Theorem 1 [Equivalent definitions of CMCPV]. 

For the PMNS mixing matrix U, the following statements are equivalent: 

(a) , it has the CMCPV (cf. Definition 1 ); 

(b) . for any parametrization of U, the general condition 01 holds; 

(c) . in the standard parametrization, (|d|, ft, 3 ) = holds; 

(d) . it has the form of RCE (up to rephasing); 

(e) . each column triangle reaches the maximal area with its e-side fixed; 

(f) . each column triangle is an isosceles triangle. 


After setting up the above preliminaries, we are ready to study the origin of CMCPV in flavor symmetry in the 
next section. 


3. Origin of Constrained Maximal CP Violation 

In this Section, we will trace CMCPV to the “residual symmetries”, i.e., the subgroups of the original flavor 
symmetry group that remain intact after the full group is broken. 

Consider that the flavor symmetry group G is broken down to two residual symmetries G v and G ( for neutrinos 
and charged leptons, respectively. They are debited as follows, 

[g : {S\S t M v S =M v }; 

G -> i (10) 

(G f : {T\T'M ( M]T = 

where M v is the Majorana mass matrix of neutrinos, and M f M] is the effective mass matrix of left-handed charged 
leptons. Thus, the mixing matrices U Y and U f (which diagonalize M v and M f M\ , respectively) are directly deter¬ 
mined by S and T lt2()h . 

ulsu v = D s , 

U]TU C = D t . 

4 


( 11 ) 




Here the matrices D s and D T are diagonal matrices. Since Eq. (11 lb demonstrates a direct connection between flavor 


symmetry and lepton mixings, it also attracted extensive studies 112 ll 221 23|,|24] via the approach of symmetry and 
group theory, without resorting to explicit mass matrices or a fundamental Lagrangian. With this general mass- 
independent approach, we will analyze the origin of CMCPV. 

Roughly speaking, our theorem states that the CMCPV can be realized if the residual symmetries are real. In 
rigorous manner, we formulate this theorem in the following form. 


Theorem 2 [Origin of CMCPV]. 

If the residual flavor symmetries in the lepton sector (including charged leptons and neutrinos) are real and fully 
determine the mixing pattern, then the CMCPV always holds, up to a few minor exceptions: 

(i) . one of the three mixing angles in the PMNS matrix is zero; 

(ii) . neutrinos are not Majorana fermions; 

(iii) . the residual symmetry for charged leptons is a Klein group, i.e., G [ — Z 1 ®Z 1 . 


It is clear that the exception (i) is already excluded by current oscillation data, and the exception fii) is not a concern 
for most neutrino theories. The exception (iii) is less trivial, but can be easily evaded in model-buildings. Besides, 
in Theorem 2, for the residual symmetries being real, we mean that there always exists a basis under which these 
symmetries become real. 

To illustrate this theorem explicitly, we first consider a simple (unrealistic) example. A rotation of 120° around 
the axis (1,1, l) r can be represented by 


R = 


0 10 ' 
0 0 1 

1 0 0 , 


( 12 ) 


Suppose that this R corresponds to the residual symmetry group G f . Thus, we have 


U ] W RU W = diag(l, to, to 2 ), 


(13) 


where U w is the Wolfenstein mixing matrix defined in Eq. (QJ. According to the relation (11 11 1, we have U ( = U w , 
which shows that Uj exhibits RCE. If we further assume that the neutrino mass matrix is diagonal (and thus U v = I), 
then the PMNS matrix U = t/J U v gives, ([d[, ft, 3 ) = (|, |), because Theorem 1 states that the RCE always leads to 
CMCPV. 

Two remarks are in order for this example. One is that both residual symmetries are real. This is explicit for G ( . 
The form of G v is, G v D {diag(l, 1,-1), diag(l, -1,1)}, since neutrino mass matrix is taken to be diagonal. The other 
point is that t/J exhibits RCE. 

In general, the validity of Theorem 2 (excluding its exceptions) implies the following important points: 


1). real G v leads to real U v , which will be explicitly proven in Appendix B.2[ 


2). real G, leads to RCE in t/J, which will be explicitly proven in Appendix B.3 


3). if t/J has RCE and U v is real, then the PMNS matrix U = U\U Y has RCE. 


The three points above are combined to prove Theorem2. The last point is just based on the Lemma given above 
Eq. and the first two points can be understood by the following reasoning. (We delegate the mathematical proofs 
to Appendix B ) Since both G ( and G v contain only real transformations, they can be geometrically regarded as 
rotations in 3-dimensional Euclidean space. (Here a trivial minus sign between the determinants of SO(3) and 0(3) 
does not matter.) For such a rotation represented by a matrix R, the rotation axis is one of its eigenvectors with 
the corresponding eigenvalue equal to 1. The remaining two eigenvectors must be complex conjugate to each other, 
which is a general property of SO(3) matrices. (The explicit forms of the two eigenvectors are given in Appendix 
B.3.) Hence, if R e G ( , then the eigenvectors of R are the columns of U f , which implies two columns in U f are 


conjugate to each other, and thus t/J has RCE. 

There are differences in the neutrino sector, because we consider the neutrinos as Majorana particles here. Hence, 
the residual symmetry has to be constructed with Zf s, i.e., G v = Z 2 ® Z 2 , which geometrically correspond to two 
rotations of 180°. These are special rotations in the sense that only such rotations may commute with rotations around 
different axes. For rotations of 180°, the eigenvalues are (1, -1, -1), cf. Eq. ( IB. 7k Due to a partial degeneracy of 
the eigenvalues, the neutrino mass matrix M v should be determined by two Z 0 -rotations with orthogonal axes. Each 
axis determines one column of U v , so U v only contains real column vectors. This in turn implies that G f cannot be 
Z 2 ® Z 2 , which is the exceptional case (iii) pointed out in Theorem2: if G f = Z 2 ® Z 2 , then U ( will be real and no 


5 
















CP violation exists. Now, it is also easy to understand why Theorem2 requires neutrinos to be Majorana fermions, 
since the symmetries Z 2 ® Z, are needed for G v . 

Theorem 2 further leads to a series of corollaries which we will discuss as follows. 

Corollary A [0(3) Subgroups]: 

If an 0(3) subgroup G contains sufficient residual symmetries that can fully determine a mixing matrix, then it leads 
to the CMCPV after avoiding the three exceptions listed in Theorem2. 

This is manifest because the constraint which requires the residual symmetries to be subgroups of 0(3) makes G ( 
and G v automatically real. According to Theorem 2, this leads to the CMCPV. Examples of such residual symmetries 
include popular groups like A 4 , S 4 , and ,4_, corresponding to tetrahedral, octahedral, and icosahedral symmetries, 
respectively. 

We should comment on the phrase “sufficient residual symmetries” in Corollary A. As is well-known, the maximal 
residual symmetries in the charged lepton and neutrino sectors are U(l) ® U(l) ® U( 1) and Z, ® Z, ® Z,, respec¬ 
tively |[20h . But when seeking flavor groups to unify the residual symmetries, it is not necessary to cover those large 
groups. For charged leptons, the minimal choice is to take a Z 3 subgroup from those U{ 1 )’s, which is in fact sufficient 
to determine the mixing U ( . For the neutrino sector, the minimal sufficient residual symmetry should be Z, ® Z,. 
So, this smaller set of residual symmetries should be included in the 0(3) subgroup for the Corollary A. 

However, in some models, especially those based on A 4 , sometimes the flavor group does not contain sufficient 
residual symmetries, and the so-called accidental symmetries are present to fully determine the mixings. Those ac¬ 
cidental symmetries depend on the detailed dynamics of the model (instead of the flavor group), so the Corollary A 
does not apply. But, if the accidental symmetry is a real symmetry, then Theorem 2 applies and there is still CMCPV. 


Corollary B [ Real M v ]: 

If G ( is real and M v is real or M v can be written as M v = z 1 1 + Z 2 M V , where M v and I are real and identity 
matrices, respectively, and (zj, z 2 ) tire complex numbers, then there is CMCPV, after evading the three exceptions 
listed in Theorem 2. 


We first consider the case that M v is real. Then, as a real symmetric matrix, M v can be diagonalized by a real 
orthogonal matrix, which implies U v and G v are real. Hence, according to Theorem 2, we have CMCPV. Multiplying 
M v by an overall complex phase will not change U v . Thus, if M v = z n M v with real M v , then this means that M v 
is essentially real, up to an overall complex phase factor. Hence, this case also leads to CMCPV. Next, consider 
M v = z\I + M v with real M v . This means that M v is real up to subtracting a constant from all diagonal elements. In 
this case, for S (e G v ) satisfying S T M V S = M v , we have S T M V S = S T (z x I + M V )S = z x I + M v , which shows that 
M v is also invariant under S. Hence, M v has invariance under real G v and there is CMCPV. Finally, combining the 
two cases above, we have thus proven the Corollary B for the general form M v = z x I + Z 0 M V . 

The form M v = z, { I + z 0 M v has important applications in model buildings for CMCPV. Typically, for building 
flavor symmetry models, at least one flavon <!> is introduced to couple with neutrinos v and forms a Yukawa term 
vv<p , which contributes to neutrino masses if the vacuum expectation values (VEV) (<p) + 0, where f is a scalar 
field acting as a multiplet under the flavor symmetry. Neutrinos are commonly considered as flavor triplets in many 
models, so a vv term or vvf term will usually show up, where f is a flavor singlet. These terms will contribute to 
M v as a diagonal mass term z x I , where Zj is complex because the coefficients (Yukawa couplings) of these terms 
are complex in general. The vvf term will contribute in a form of z~,M v if ((b) is real (up to an overall complex 
phase) and the Clebsch-Gordan (CG) coefficients for the term are real. 

Note that here we only consider the case with all 3 generations of neutrinos unified into a triplet of the flavor 
group. Otherwise, the vv term would not be diagonal^ The flavon <b can be in any non-trivial representation. We 
should point out that both real (<p) and real CG coefficients are very common in many groups. For instance, in the 
3-dimensional representation of A 4 , the CG coefficients for 3 ® 3 ® 3 —» 1 are real in both the real basis (used for 
example in 12511 and the complex basis (used for example in [26]). This corollary does not apply to groups with 
inherent complex CG coefficients, like T' [27] or A(27) [28]. As for real {(b) , if <b is a real scalar field by definition, 
then (<f>) is real. If 4> has to be a complex field, then as known from minimization of scalar potentials in many models, 
it is still common to have VEV alignment according to (<p) oc (1,1,1), (1,0,0), etc., which is real. If (<p) is however 
complex, then in general CMCPV does not follow. For those “real vv<p models”, where only the vv(b term makes a 


Tor certain groups, such as A(3 n 2 ) or E(3n 3 ) with n ^ 3 , two triplets cannot form a singlet. Those are certain subgroups of SU(3), to be 
precise, subgroups with faithful irreducible 3-dimensional representation whose determinant equals 1, that have complex representations and are 
not subgroups of SO(3). In this case, the vv term is absent, which means that in the general form M v = z l I + Z 2 M V only z 2 M v exists. So the 
problem becomes simpler. 
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non-trivial contribution (not proportional to the unit matrix) to M v , we have the following corollary. 

Corollary C (Real vvcf> Models ]: 

Real wcp models always lead to CMCPV, if the three minor exceptions listed in Theorem 2 do not happen. 

For a demonstration of the above general discussion, we will build a simple “real vv<p model” in the following 
Section l4~3l 


4. Applications 

In this section, we apply our theorems and corollaries to various situations and understand why CMCPV is realized 
in certain cases. We will illustrate how to achieve the CMCPV in model buildings. There are extensive recent literature 
| 29h studying specific models of (|5|, 0 23 ) = (|, |). 


4.1. p — t Reflection Symmetry 

The p-r reflection symmetry was studied before QSIIQI , which is sometimes also called the generalized p-r 
symmetry. This symmetry is defined in the flavor eigenbasis (with M f diagonal), 
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e 




(14) 


Imposing this symmetry leads to the following form of the neutrino and lepton mass matrices. 
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(15) 


where M v is symmetric and Mj = M ( M', is diagonal. Note that the elements r l , are real, but z 12 are complex 
in general. The operation dT4b will transform v r M v v to its Hermitian conjugate. We can directly check that the 
Lagrangian term, _£ D v T M v v + h.c., is invariant under the transformation (fl4l > if and only if M v takes the form of 
Eq. O. 

Let us make a transformation, 

M v = U ( M y U t ( , M) = U ( M]U], (16) 


with 


Thus, we derive 
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(18) 


where we have defined notations, z.j = z ;1 + tip , (j = 1,2), and a = m ~, b ± = \{m 2 ^ ± m 2 ) . The quantities (Zj V Zj 2 ) 
and (a, b ± ) are all real. Note that M v is a real matrix, and the charged lepton sector has an SO(2) residual symmetry. 


R = 


10 O' 
0 cos 8 sin 8 
0 - sin 8 cos 8 t 


(19) 


This is because RM^R^ = Mj holds for 8 e [0, lit). Since M v and G f are all real, this will lead to CMCPV 
according to our Corollary B. Note that real M v implies that G v is real. 

The p-r reflection symmetry is certainly not the only possibility to generate CMCPV. From Eq. (IT8l> . we see that 
M v is real, while Corollary B shows that M v can have a more general form including complex numbers. Hence, the 
p-r reflection symmetry is just a special case of real residual symmetries, although this is not manifest before the 
transformation of basis in Eq. ( 1 1 6k 
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4.2. CMCPV from Geometrical Symmetry Breaking 


As another example illustrating our theorem, we revisit a model from Ref. |23|], which predicted <5| = ? and 
0 23 = | (as well as 0 |3 - | - 0 p ). We will show that this model fulfills the criteria for CMCPV. 

This model identifies a Z 4 rotation around the jr-axis as G f , and the product reflections Z, ®Z 0 as G,„ where one 
Z 2 reflects y —> —y and the other Z 2 transforms (x, z) —> -(z, x). These rotations are subgroups of the octahedral 
symmetry O h , and can be shown by simple geometrical picture. This group setting generates the bimaximal mixing, 
0 |2 =023 = f and 0 |3 = 0. The necessary deviation from this leading order scheme was generated by slightly 
tilting the axis of Z 4 rotation by a small angle (defined as y2e ) that turns out to be related to nonzero 0 |3 . We also 
verified that this geometrical symmetry breaking can arise from certain flavon models. For example, we may set up a 
concrete realization, where a flavor triplet (p is responsible for mass-generation of the charged leptons and the Yukawa 
terms involving if) are SO(3) symmetric in the 3-dimensional flavor space. With these, the geometrical breaking is 
connected to the VEV misalignment of flavons. After the axis tilt, the residual symmetry of charged lepton mass 
matrix is represented by 12311 , 

1 -2e 0 


R ( = 


0 

2e 


-1 

0 


+ 0(e 2 ). 


( 20 ) 


The neutrino mass matrix is still invariant under the original reflections Z,, ® Z 2 , as represented by 



' 0 

0 

-1 ' 


' i 

0 

O' 

*vl - 

0 

1 

0 


0 

-1 

0 


l-l 

0 

0 J 


lo 

0 

1J 


( 21 ) 


Since all the residual symmetries are real. Theorem 2 applies and the model should realize CMCPV. This is indeed 
the case, as found in Ref. | 


4.3. A Real vv(p Model 

As stated in Corollary C, the real vv(f> models should always produce CMCPV. In the following, we will build 
such a model as an explicit illustration. 

We use A 4 <g> Z 2 as flavor symmetry group and introduce two scalar fields <// and <p v , in addition to the SM 
Higgs doublet H. The relevant particle content of this model is summarized in Table[T| The Lagrangian for the 
lepton-neutrino sector contains, 

-C D y e (Ltp e )He c + y^L^YH/Z + y T (L<p { )" Ht c 

+ y vl (LL)HH + y vl {<p v LL)HH + h.c., (22) 

where L stands for the left-handed lepton doublet of SU(2) L and H is the Higgs doublet. Since <f> c and <// do 
not have any charge other than the Z 0 assignment, they can be real fields. Consider that they acquire the following 
alignment of VEVs, 

<0 v >oc(l, e 2 , e 3 ), <0 f )oc (1,1,1), (23) 

where similar to we have introduced a small perturbation on the usual VEV alignment in (<p v ), with e 2 , e 3 c 1. 
Similar VEV alignment was already considered in the literature a , but its further elaboration is irrelevant to the 
current illustration purpose of realizing CMCPV on the ground of residual symmetries li20ll ; it is also fully beyond the 
main goal of this short Letter. For e 2 = e 3 = 0, one would obtain the tri-bimaximal mixing; and the small (e 2 , e 3 ) 
should produce the necessary corrections a. Note that for e 2 - e 3 = 0, although {(p v ) is real as required by 
our definition of real vvcp models, one of the mixing angles, 0 13 , is zero, which just matches the exception-(i) of our 
Theorem2. Hence, CMCPV does not follow in this case. For e 23 4 0. the charged leptons and neutrinos acquire 
masses as follows, 



' a 

b 

c 


' d 

e 3 

f 2 ' 

Mf oc 

a 

bm 2 

C(j) 

, M v OC 

e 3 

d 

1 



bo 

cur , 


- e 2 

1 

d, 


where the mass parameters (a, b, c, d) are complex in general. This type of lepton and neutrino mass matrices 
are often studied in the literature im For instance, diagonalizing the lepton mass matrix M' ( M [ gives the mass- 
eigenvalues (m e , m fI , m T ) oc (| a\, \b\, |c|). [This also shows that using the observed mass values (m e , m , m T ) does not 
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Table 1: Particle content of the A 4 <S> Z 2 model. 


Groups 

L 

(e c , n c , t c ) 

/ 

r 

H 


3 

(1, 1", 1') 

3 

3 

1 

z 2 

-1 

1 

-1 

1 

1 

SU(2), 

2 

1 

1 

1 

2 


fully fix the parameters (a, b, c ) themselves; while inputting the model parameters (a, b, c) can fully accommodate 
the observed lepton masses.] The focus of our paper is on the origin of leptonic mixings from flavor symmetry. Thus, 
by diagonalizing M ( and M v , we derive the following lepton and neutrino mixing matrices, 


and 


( f 2 +e 3 

1 

-e 2 +e 3 3 

V3 

V2 

V2 

~ 6 2 

V2 

1 


1 

( V2 

f 3 

V5 / 


if 11 1 ' 

U r = —— 1 a) 2 u> 

^ , 1 oj or, 


(25) 


(26) 


From U = UlU Y , it is straightforward to extract the PMNS parameters in the standard parametrization, 


n 

23 4 ’ 


n 
2 ’ 


(27a) 


sin# 13 ^ 



tan0 p =* 


V 2 (l-e 2 -e 3 ) 

2 + e 2 + e 3 


(27b) 


These results show that, apart from model-specific deviations of (0 13 , 0 p ) from their tri-bimaximal values, we have 
realized the CMCPV, as expected from Corollary C. From Eq. (I27bl) . we can determine the perturbative parameters 
(e 2 , e 3 ) in terms of (sin$ 13 , tan0 p ) via 





1 - V2 tanfip 

2 +V2 tanffp 


e 3 



+ 


1 - V2 tan0 12 

2 + V2 tan0j 2 


(28) 


Taking 0 |3 =* 9° and 9 V - 34°, we derive (e 2 , e 3 ) - (-0.18, 0.21). 

In general, we can extend the real vvcp models to type-I neutrino seesaw. In this case, we may introduce three 
right-handed neutrinos v R in the 3-dimensional representation of A 4 . Thus, the neutrino Dirac mass matrix will be 
proportional to unit matrix, m D oc I , while the heavy Majorana mass matrix M R shares similar structure with M v 
in Eq. (l24l >. Hence, we find that the seesaw mass matrix of light neutrinos M v oc M R 1 . 


5. Conclusions 

In this work, we stressed that a general parametrization-independent definition of the maximal CP violation should 
be constructed in terms of Jarlskog invariant J , rather than the CP angle 6 (which is rephasing non-invariant). We 
pointed out that naively maximizing J without constraint is already excluded by oscillation data. We further demon¬ 
strated the crucial importance of introducing the new concept of constrained maximal CP violation (CMCPV) for 
studying the viable maximal CP violation. For this purpose, we constructed CMCPV in the Definition 1, and formu¬ 
lated it by a set of equivalent ways, as summarized in Theorem 1 (Sec. 2). We derived the parametrization-independent 
realization of the CMCPV via solution (0), which was proven to be the maximum of Jarlskog invariant under a min¬ 
imal constraint on the PMNS matrix U (Sec. 2 and Appendix A). We found that the CMCPV just corresponds to 
(|5[, 4,,) = (f, |) in the standard parametrization of the PMNS matrix (0. In Sec. 3 and AppendixB, we proved 
Theorem 2, stating that if the residual symmetries of neutrinos and charged leptons are real, then the CMCPV should 
be realized, up to a few minor exceptions. It was shown that the conditions for CMCPV are actually quite common, 
and we presented several sample models in Sec. 4, demonstrating that in particular the often considered /j-t reflection 
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symmetry is a special case of our theorem. We also note that the current formulation cannot be naively applied to the 
quark sector. The reason is that our Theorem 1 proves RCE to be essential for the CMCPV, but RCE cannot hold for 
the CKM mixing matrix due to experimental data. Namely, any two rows (or columns) in the CKM matrix cannot be 
conjugate to each other (up to rephasing). 

If indeed the values of 8 =* — | and 0 23 — f continue to be favored by neutrino data, our general theorems 
and corollaries of CMCPV should be important, and provide strong guidelines for the model buildings with flavor 
symmetry. 


Appendix A. Parametrization-Independent Solution of CMCPV 

In this Appendix, we derive the general solution of CMCPV by using the manifestly parametrization-independent 
formula of Jarlskog invariant £j}. 

Following our Definition 1 for CMCPV, we can use Eq. (I3bb to derive the extremal conditions of Jarlskog invariant 
respect to \U M i\ and \U^\ by fixing \U e \\ and \U e -}\. Thus, we have 


dJ 2 

1 



dz 

= xyw - - (1 - y) [(1 - y) z + (1 - x) w - (1 - x - y)J 

= 0 , 

(A. la) 

dJ 2 

1 



dw 

= xyz- -(l-x)[(l-y)z + (l-x)w-(l-x-y)J 

= 0 , 

(A. lb) 


where for convenience we have used the notations, (x, y, z, w) = (|f/ e i| 2 , \U e i\ 2 , |t/^i| 2 , It/^l 2 ). From the extremal 
conditions (IA. laMA.lbb , we deduce the solutions, 

z=^(l-*), vv=i(l-y). (A.2) 


Hence, we have 


Kj\ 2 = {{'-Wef)’ O' =1,2,3), (A.3) 

3 

where we have used the unitarity condition for the second row, ^ |(/ w j 2 = 1. With Eq. (IA.3I) and making use of the 

j= i 

unitarity conditions for each column of the mixing matrix U , we further deduce 

\U Tj \ 2 = 1 ( 1 - 1 ^/), 0 =1,2,3). (A.4) 

Finally, comparing Eqs. (IA.3b abd (IA.4b . we arrive at 

\uj = I U Tj \ 2 = |(l-| u ej \ 2 ), 0= 1,2,3). (A.5) 

This just reproduces the Eq. ©, which we presented in the text. 

Next, we prove that the above extremal solution (IA.2b or (lA.3b indeed corresponds to a maximum of Jarlskog 
invariant. For this purpose, we compute the second derivatives of the squared Jarlskog invariant respect to (z, w) , 

c^ 2 )" = -|a-y) 2 , v 2 y;, w =-\a-x)\ u 2 )" v = uX z = -\v-x-y-xy). (a.6> 

Then, we inspect the eigenvalues of the 2x2 matrix {(/ 2 )"}, whose elements are given by Eq. (IA.6I >. The eigenvalues 

{/Ij, /1 2 ) satisfy the following quadratic eigenvalue equation. 


A 2 - BA + C = 0, 

(A.7a) 

B = (J 2 )" z + (J 2 )'ww ~ -^[d-x) 2 +(l-y) 2 ] < 0, 

(A. 7b) 

C - (/ 2 )"(/ 2 )" w -[(/ 2 )"J 2 = xy(l-x-y) >0, 

(A.7c) 
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where we have, 1-x-y- 1- \U e \\ 2 - |t/ C 3 | 2 = [ t/ e 2 1 2 >0, due to the unitarity condition on the first row. Thus, we 
have the two eigenvalues obey A 1 + A 2 = B < 0 and A l A 1 = C > 0. This means that the two eigenvalues of {(/ 2 )"} 
are both negative, A v A-, < 0. Hence, we conclude that the extremal solution (IA.2b or (IA.3b is indeed the maximum 
of the Jarlskog invariant (under the constraint on the first row of U ), and provides the parametrization-independent 
realization of the CMCPV as given in our Definition 1. 

Finally, using the parametrization-independent solution (lA.5b or ([4} of CMCPV, we can readily derive the ex¬ 
plicit realization of CMCPV under the standard parametrization (0. From the first equality of Eq. (1A.5b . we have 
two independent conditions \U M \\ = \U T \\ and \U^\ = \U T j\, which take the following forms under the standard 
parametrization 0. 

| i 12 c 23 + c 12' s 23' s 13 el<5 | = C 12 c 23 S 13 eI<5 | > 

| S 23 C 13| = | C 23 C 131 ' 

The condition (IA.8bl > leads to s 23 = c 23 and thus 0 O3 = |. Given this, we can rewrite (IA.8al > as 

kl2 +c 12*13 eI 1 = k2- c 12 s 13 e ‘1- ( A -9) 

Since c 12 s 13 4 0, this must require cos <5 - 0, i.e., |<5| = | . Hence, the explicit realization of our CMCPV under the 
standard parametrization (0 just gives (|d|, d 23 ) = (|, |), which we mentioned in the text above Eq. 0. 

Appendix B. Proofs 

In this Appendix, we present proofs that are needed to establish the Lemma given after Eq. 0 and the main 
Theorem2 given in Sec. [3] 

Appendix B.l. RCE is Invariant under Right-handed Real Transformations 

For a unitary matrix V with the form of “row conjugation equality” (RCE) and a real orthogonal matrix R , we 
need to prove that V' = VR still has RCE. The proof is straightforward. Defining the elements of these matrices, 

V = (u u ), R = (r u ), V'=(u' u ), (B.l) 

we have 

«y =Z M *V (R2) 

k 

Note that the matrix elements (u lk ) and (r k] ) (with k, j - 1,2,3) are real numbers from the start. The RCE feature 
of matrix V gives, (u 2j )* = « 3j for j - 1,2,3. This implies 

u[j = ^ u \k r kj = rea l numbers, (B.3) 

k 


(A.8a) 
(A.8b) 


and 


(U 2?* = X “2 k r kj = 2 U3kn ’ = "3 j ■ •'B' 4 ) 

k k 

We have thus proven explicitly that RCE is invariant under right-handed real transformations. 


Appendix B.2. Real G v Leads to Real U v 

Consider Majorana neutrinos with residual symmetry G v = Z 2 ® Z 2 . In the following, we will prove that a real 
G v leads to real U v , and vice versa. 

Let us set S to be a 3x3 unitary matrix, which is real and is a Z, transformation (i.e., .S’ 2 - I ). As S is real, it 
follows that S k = S T , and the unitarity condition S S ' = I implies that the real matrix S is orthogonal, S S 7 - I. 
Without losing generality, we set S 6 SO(3). Hence, S is a rotation in 3-dimensional Euclidean space. Furthermore, 
since S 2 — I , it must be a 180°-rotation. 
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For G v - Z, (g) Z,, we may use S l and S 2 to represent the transformations of the two Z,’s, respectively. Thus, 
[5 j, 5,] = 0 should hold, which implies that their rotation axes must be orthogonal. Hence, geometrically G v 
contains two 180°-rotations with orthogonal axes. These two axes can be represented by two normalized real vectors 
Vj and v 2 with 


Si 

s 2 


S 1 1 '2 = ~ V 2 
S 2^2 — v 2 , 


(B.5) 


where Vj and v 2 are column vectors, of the 3x1 matrix form. Taking v, = VjXv 2 and U v = (v^, v 2 , v 3 ), we see that 
U v is a real matrix and can diagonalize S l and S 2 simultaneously in the way given by Eq. i flTTl . 

Therefore, if G v D {S l , S 2 } contains only real matrices, then U v is real. The converse proposition that a real U v 
leads to real G v is also true, and can be readily proven. 


Appendix B.3. Real G { Leads to Complex Uwith RCE 


We need to prove that any SO(3) matrix R can be diagonalized by U r R U r , where the unitary matrix U R contains 
one real column and two other columns which are complex conjugate to each other. This can be explicitly proven as 
below. 

The most general rotation in 3d Euclidean space which rotates the space around an axis n = (n l , n 2 , n 3 ) T by an 
angle (f> is 12411 , 


n 3 +c(n 2 +n^j 

R(n,<f>) = (1 — — s« 3 c+nt;—cn 2 sn , 

sn 2 +(l— cjUjUj — s/ij +(1 — c)n 2 « 3 ^ . ,* 3 «,#» 3 

where n • n = 1 and (s, c) = (sin (f>, cos </>). We can directly verify that this matrix is diagonalized as 


(1 — c)n 1 n 2 +sn 3 — sn,+(l— c)nj« 3 ^ 

-(1 — c)n 2 n 

ji 


T 
c+ni 


-cn 


2 3 
2 


(B.6) 


where 


u\ru r 


10 O' 
0 c+ii 0 
0 0 c-is / 


VT^ 

Vl-”I j 

”i V2 

V2 

n.n~- in. 

n l n-,+in ,, 

Z Sj2a-ri\) 

V2(l-nj) 

«l«3 +i «2 

in 7 

J V 2 d- n ?) 

y 2 (i-«j) J 


(B.7) 


(B.8) 


We see explicitly that the first column of U R is real, and the second and third columns are conjugate to each other, 
i.e., U' R has RCE. Hence, if R e G f , then u\ = U R has RCE. 
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Note added: While we were finalizing the present paper. Ref. 13211 appeared on arXiv, which has some partial overlap. 
It was also pointed out there that p-r reflection symmetry can be generated by discrete residual subgroups of 0(3). In 
Sec. 14. II we explicitly showed that with a proper basis transformation the p-r reflection symmetry is actually a real 
symmetry. Our general theorems are independent and complementary to [i32[], and we presented a set of equivalent 
formulations for the CMCPV as well as its parametrization-independent realization. 
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